Protocol for high fidehty readout in the photon blockade regime of circuit QED 
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The driven-damped Jaynes-Cummings model in the regime of strong couphng is found to exhibit 
a coexistence between the quantum photon blockaded state and a quasi- coherent bright state. We 
characterize the slow time scales and the basin of attraction of these metastable states using full 
quantum simulations. This form of bist ability can be useful for implementing a qubit readout 
scheme that does not require additional circuit elements. We propose a coherent control sequence 
that makes use of a simple linear chirp of drive amplitude and frequency as well as qubit frequency. 
By optimizing the parameters of the system and the control pulse we demonstrate theoretically 
very high readout fidelities (> 98%) and high contrast, with experimentally realistic parameters for 
qubits implemented in the circuit QED architecture. 
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Qubit readout in solid state systems is an open prob- 
lem, which is currently the subject of intensive experi- 
mental and theoretical research. High fidelity single shot 
readout is an important component for the successful im- 
plementation of quantum information protocols, such as 
measurement based error correction codes [1 as well as 
for closing the measurement loophole in Bell tests [2l |3] . 
For measurements where the observed pointer state de- 
pends linearly on the qubit state, for example dispersive 
readout in circuit QED (cQED) [4], there exists a uni- 
fied theoretical understanding [F. Experimentally, these 
schemes require a following amplifier of high gain and 
low noise, spurring the development of quantum limited 
amplifiers p]. However, the highest demonstrated fideli- 
ties to date rely on nonlinear measurement schemes with 
qubit dependent latching into a clearly distinguishable 
state, e.g. Josephson Bifurcation Amplifier (JBA) as well 
as optimized readout of phase qubits [3l[7]- For this class, 
during the measurement, the system evolves under the in- 
fluence of time varying external fields and nonlinear dy- 
namics, ultimately projecting the qubit state. The space 
for design and control parameters is very large, and the 
dependence of the readout fidelity on them is highly non- 
trivial. Therefore the optimization is difficult and does 
not posses a generic structure. 

We propose a coherent control based approach to the 
readout of a qubit that is strongly coupled to a cavity, 
based on an existing cQED architecture, but not neces- 
sarily limited to it. This approach is in the spirit of the 
latching readout schemes, but it differs in that the source 
of the nonlinearity is the Jaynes-Cummings (JC) interac- 
tion. When the qubit is brought into resonance with the 
cavity mode, the strong anharmonicity of the JC ladder 
of dressed states can prevent the excitation of the system 
even in the presence of a strong drive, a quantum phe- 
nomenon known as photon blockade [8 . However, due to 
fact that the JC anharmonicity is diminishing with the 
excitation number, we find a form of bist ability, where 
highly excited quasi-coherent states (QCS) co-exist with 
the blockaded dim states (Fig. [T]). In order to make use 
of this bist ability to read out the qubit, it is necessary 
to solve the coherent control problem of selective popu- 



lation transfer, which is how to steer the system towards 
either the dim state or the QCS, depending on the ini- 
tial state of the qubit (Fig. |3|. This selective dynamical 
mapping of the qubit state to the dim/bright states con- 
stitutes the readout scheme. It is potentially of high con- 
trast, and hence robust against external amplifier imper- 
fections. An advantage of this readout is that it uses no 
additional components beyond the qubit and the cavity, 
the latter already present as part of the cQED architec- 
ture. Based on a full quantum simulation which includes 
dissipation of the qubit and the cavity (we ignore pure 
dephasing [9 ) , we predict that implementing this scheme 
should yield very high fidelities between 90% and 98% for 
a typical range of realistic cQED parameters (Fig. [4| . 
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Figure 1. A schematic of the dim state (left histogram) in the 
quantum part of the Jaynes-Cummings (+) manifold and the 
QCS in the semi-classical part (right histogram). Application 
of a readout pulse sequence dynamically maps the initial qubit 
state into one of the states with high fidelity. The inset shows 
the lifetimes of several QCS with increasing drive strength 
and mean photon number, while for these parameters the dim 
state is photon blockaded (the coexistence regime), and the 
qubit decay (Ti) has a distinct influence on the lifetime (rb) 
of the QCS. 
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We consider the Jaynes-Cummings model with drive 
and dissipation (h = 1) 

H = cjca'''a+ Ycr^+^(cr+a+aV_)+^(t)(a+a''')+i^^+i7;^ 

(1) 

where ^(t) is the time-dependent drive of the cavity, g 
is the cavity-qubit couphng, and H^ ,^ represent the cou- 
phng to the qubit and cavity baths, respectively. When 
the system is initialized in the ground state, there is a 
range of drive strengths for which the system will remain 
blockaded from excitations out of the ground state. How- 
ever, since the anharmonicity of the JC ladder decreases 
with excitation number, the transition frequency for ex- 
citations between adjacent levels ultimately approaches 
the bare cavity frequency. Qualitatively, when the ex- 
citation level n is such that the anharmonicity becomes 
smaller than the linewidth a^, we expect the state dy- 
namics to be semiclassical, similar to a driven-damped 
harmonic oscillator \10^. More specifically, in order to 
support a coherent wave packet centered around level n, 
with a standard deviation of ^/n^ the difference of tran- 
sition frequencies across the wavepacket has to be of the 
order of the linewidth Hi. This approximate criterion for 
a minimal n can be written as 0Jn-\-2a — ^n-2a ^ ^ where 
^n±2a are the ladder transition frequencies, positioned 
2(7 = 2^/n above and below the mean level n. Quanti- 
tatively, we find in simulations that the lifetime of the 
QCS is long but finite, and increases with the amplitude 
of the drive (see insert in Fig. [T]). As we explain below, 
we find that low lying QCS states (n = 20) are the most 
effective for optimizing the overall readout fidelity. Note 
that the JC ladder consists of two manifolds (originating 
from the degeneracy of the bare states |^, n+1) and |e, n)) 
denoted by (±), and we will always refer to a states oc- 
cupying one manifold since the drive is off-resonant with 
respect to the other manifold. Transitions between man- 
ifolds contribute to the decay of the QCS to the dim 
state. Such transitions can be induced by the drive but 
their rate is smaller by a factor of 0(n~^/^) compared to 
the rate of transitions inside the same manifold. An addi- 
tional source of inter-manifold transitions are decay (Ti ) 
and pure qubit dephasing (T^^), whose effects in the pres- 
ence of drive were studied in the context of the dispersive 
regime [TT]. Indeed, as we see in Fig. [l] changing Ti has 
a noticeable effect on the QCS lifetimes. For very large 
n, these processes become ineffective for inducing decay 
of QCS, since then the difference between manifolds exci- 
tation frequencies becomes smaller than and therefore 
the drive effectively drives both manifolds. For supercon- 
ducting transmon qubits Ti is the dominant decay pro- 
cess, and we show its effect on the overall fidelity in Fig.|4] 
The QCS exist with drive amplitudes where the ground 
state is photon blockaded, giving rise to a dynamical 
bistability between quantum and semi-classical parts of 
the JC ladder. Indeed, we see that there is a basin of 
attraction for states initialized as coherent wave pack- 
ets to persist as QCS, and we characterize it according 
to the probability of the state to decay on the timescale 
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Figure 2. Probability for the QCS to decay after being ini- 
tialized as coherent state wave packet |ae*^) and driven for 
a time . The equal probability contours trace out two 
basins of attraction: states initialized inside the low decay 
probability contour {Pd < 1%) end up long lived (n ^ 
whereas states initialized left of Pd = 90% quickly decay to 
the ground state and remain photon blockaded. The pa- 
rameters are g/27r = lOOMHz, k/27v = 4.05MHz, ^/27v = 
9.9MHz,(cjd - uJc)/27T = 12.3MHz, (cjc - cjg)/27r = 9.7MHz, 
Ti = 1591ns {(jJd is the drive frequency). 



Tb':^ K,~^. In Fig. [2] we plot the contours of equal proba- 
bility of the QCS to decay to a manifold of states close to 
the ground states, after a time given that it was ini- 
tialized with a certain amplitude {a) and phase {6). We 
see a large region supporting QCS, and the phase sen- 
sitivity can be understood qualitatively from the time 
dependent simulations: a mismatch between the phases 
of the drive and initial coherent state causes ringing of 
the wave packet outside of the basin into the too anhar- 
monic part of the ladder from which it cannot recover. In 
addition to the existence of this basin, the anharmonicity 
acts together with the cavity decay to induce mixing of 
the QCS: even for bright states (n > 40) we extracted a 
relatively low purity of Tr(p^) < 0.5. 

In the regime of coexistence the dim quantum state 
and QCS present us with the possibility of implementing 
a high contrast readout scheme. This requires the solu- 
tion of the coherent control problem of steering the logical 
It) state to some point on within the basin of attraction 
(Fig. [2|, while keeping the \D far from the basin, in the 
manifold of dim states. In the presence of dissipation, 
the latter would quickly decay to the ground state, and 
remain there even in the presence of driving, due to the 
photon-blockade, whereas the QCS would persist for a 
long time and emit approximately n{n)i)Ti) photons. 
The standard coherent control problem of population 
transfer [12], which was also discussed recently for su- 
perconducting qubits [13] , is to maximize the probability 
Pi^f of steering the state \i) to the state |/). However, 
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Figure 3. Readout control pulse (a) Time trace of the drive 
amplitude: a fast and selective initial chirp (10ns) can se- 
lectively steer the initial state, while the qubit is detuned 
from the cavity {{ujq — ujc)/'^7r ^ 2g). It is followed by a 
slow displacement to increase contrast and lifetime of the 
latching state, while the qubit is resonant with the cavity 
{k/2ti = 2.5MHz). The drive amplitude ramp is limited so 
that the photon blockade is not broken, but the contrast is 
enhanced by additional driving at the highest drive ampli- 
tude, (b) A diagram of transition frequencies shows how the 
drive frequency chirps through the JC ladder frequencies of 
the (+) manifold, and how the manifold changes due to the 
time dependent qubit frequency, (c) Wave packet snapshots 
at selected times (indicated by bullet points on panel (b)) of 
the chirping drive frequency of panel (b) conditioned on the 
initial state of the qubit. 



here the goal is to bring the probability for selective steer- 
ing Pi^f + Pi'^f < 2 close to its theoretical maximum, 
which is an essentially different coherent control problem. 
For systems with very large anharmonicities, for example 
atomic systems it is possible to effectively implement a 
population transfer via adiabatic control schemes such as 
STIRAP [12 . The JC ladder anharmonicity is relatively 
small compared to atomic systems, and so these schemes 
are inapplicable here. 

The control pulse sequence we apply is depicted in 
Fig. |3j and consists of three parts: (1) a strong chirped 
pulse {t < tc) drives the cavity, with the qubit being de- 
tuned such that the cavity frequency is state dependent. 
This pulse maps the |^) and |t) to the dim and bright 
state basins, respectively. (2) a weaker pulse that trans- 
fers the initially created bright state to even brighter and 
longer lived states {tc < t < th)^ and (3) steady driv- 
ing for additional contrast {th < t < tf). In designing 
such a pulse sequence we have the following physical con- 
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Figure 4. Cumulative probability distributions for the num- 
ber of photons (N) emitted from the cavity during the driving 
time tf, for different qubit decay times (Ti ) , including a very 
long Ti = 15/is indicating that Ti is not limiting the readout 
fidelity. For low detection thresholds {Nth ~ 20) for distin- 
guishing It) (A^ > Nth) from \i) {N < Nth) the fidelity can 
be very high (> 98%) for realistic values of qubit decay in 
cQED (a few microseconds), and of high contrast for more 
moderate fidelities (> 90%). The distributions also show a 
almost no false positives for higher thresholds Nth > 20 (here 
th = 194ns and other parameters as for Fig. p|. 



siderations: (a) the initial fast selective chirp has to be 
optimally matched to the level structure so that the pop- 
ulation transfer and selectivity would be extremely high 
(b) it is necessary to chirp up quickly before decay pro- 
cesses become effective and result in false negative counts 
{tcK, ~ 0.16) (c) for t > tc it is necessary to drastically re- 
duce drive strength, since it reaches drive strengths which 
would break the photon blockade through multiphoton 
processes if it persisted. The piecewise linear chirp se- 
quence is fed into a full quantum simulation that includes 
decay, and the 13 parameters of the system and drive are 
optimized with respect to the total readout fidelity. 

The cumulative probability distributions to emit A^ 
photons conditioned on starting in two initial qubit states 
are plotted in Fig. [4j These distributions were optimized 
for Ti = 1/JjS and then regenerated after varying Ti, in or- 
der to depict the effect of qubit relaxation on the readout. 
There are two figures of merits from figure [4j one is that 
there exist very high fidelities T =1- P(t | |) - P{i \ t) 
exceeding 98% for a low threshold around Nf^ = 20 
even for relatively short lived qubits (Ti ^ 500ns). In 
order to take advantage of these fidelities a very low 
noise amplifier would be needed. In addition, we find 
high contrast and high fidelities (> 90%) for long lived 
qubits Ti > 1.5/is with thresholds around Nth = 150. 
These fidelities will not be limited by the noise added 
by a cryogenic HEMT amplifier with noise temperature 
of T/v ^ 5i^. For T/v > the effect of the noise can 
be mitigated by increasing tf (up to time r^). Note also 
that the limit of obtainable fidelities with this control 
scheme is not due to finite qubit lifetime, as we see from 
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the curve that was simulated for Ti = 15 /us. Another 
useful feature of these distributions is the very low level 
of false positives (red curve for qubit state |^)), for a wide 
range of thresholds, originating from the effectiveness of 
the photon blockade. 

For experimental applications it is also important to 
know how robust the fidelity is against deviations of the 
control pulse parameters from their optimal values. We 
therefore varied the parameters of the initial chirp pulse 
Sd,c = ^d- ^c, '^d,c, ^d,q =^d- ^q, ^d,q, ^, ^c, and K inde- 
pendently around their optimal values and checked the 
range for which the fidelity is above 98%. The fidelity is 
most sensitive to variations of the duration of the chirp 
pulse tc which yields a tolerance of ±10% and higher 
ranges (> ±20%) for the rest of the parameters. 

For cQED all the necessary components for the above 
scheme have been experimentally demonstrated. Strong 
qubit-cavity coupling has been demonstrated in many ex- 
periments [14-16 . Strong driving of a cavity-qubit sys- 
tem has been shown in [il7], with the system behaving 
in a predictable way, as well as photon blockade [18] and 
fast dynamical control of the qubit frequency via flux 
bias lines |T9|. In addition there is evidence both theo- 
retically and experimentally for the increasing role that 
quantum coherent control plays in the optimization of 
these systems for tasks of quantum information process- 
ing. As examples we can mention improving single qubit 
gates [20], two-qubit gates [21 , and population trans- 



fer for phase qubits [13]. We therefore believe that the 
readout scheme would be applicable for the transmon, 
although the control parameters would have to be re- 
optimized due to the effect of additional levels. 

The suggested readout scheme is different from other 
existing schemes in several aspects. Compared to disper- 
sive readout [4 it involves very nonlinear dynamics and 
could potentially exhibit much higher fidelity and con- 
trast. Even though it relies on a dynamical bistability, 
it is essentially different from the JBA and another re- 
cently proposed scheme [22], since it explicitly operates 
using the quantum photon blockade. It is important to 
stress that the above optimization of the control param- 
eters is only partial, since we have limited ourselves to 
simple linear chirps in this work. More complex modula- 
tions are certainly possible although the standard meth- 
ods for optimal control [23 may be difficult to implement 
here due to the large Hilbert space. Therefore we believe 
that an experimentally based optimization using adap- 
tive feedback control [24] might be best option, and has 
the potential to yield superior readout fidelities for higher 
detection thresholds. 
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